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Study and modeling of finite rate chemistry 
effects in turbulent non-premixed flames 

By Luc Vervisch 
1. Background and objectives 

The development of numerical models that reflect some of the most important 
features of turbulent reacting flows requires information about the behavior of key 
quantities in well defined combustion regimes. In turbulent flames, the coupling 
between turbulent and chemical processes is so strong that it is extremely difficult 
to isolate the role played by one individual physical phenomenon. Direct numerical 
simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry 
interaction in some restricted but completely defined situations. Globally, non- 
premixed flames are controlled by two limiting regimes: the fast chemistry case, 
where the turbulent flame can be pictured as a random distribution of local chemical 
equilibrium problems; and the slow chemistry case, where the chemistry integrates 
in time the turbulent fluctuations. The Damkohler number, ratio of a mechanical 
time scale to a chemical time scale, is used to distinguish between these regimes. 
Today most of the industrial computer codes are able to perform predictions in the 
hypothesis of local equilibrium chemistry using a presumed shape for the probability 
density function (hereafter pdf) of the conserved scalar (Couplan et al. 1987, Viollet 
et al. 1990). However, the finite rate chemistry situation is of great interest because 
industrial burners usually generate regimes in which, at some points, the flame is 
undergoing local extinction or at least non— equilibrium situations. Moreover, this 
variety of situations strongly influences the production of pollutants (Masri et al. 

1988). 

To quantify finite rate chemistry effect, the interaction between a non-premixed 
flame and a free decaying turbulence is studied using DNS. The attention is focused 
on the dynamic of extinction, and an attempt is made to quantify the effect of 
the reaction on the small scale mixing process (section 2). The unequal diffusivity 
effect is also addressed (section 3). Finally, a simple turbulent combustion model 
based on the DNS observations and tractable in real flow configurations is proposed 
(section 4). 

2. Finite rate chemistry effects in turbulent non-premixed flames 

Three-dimensional direct numerical simulations of non-premixed flames have been 
carried out using a high-order finite difference code (Trouve, 1991). Reacting flows 
are simulated resolving all of the scales of the turbulent motion, including effects 
due to variations in density and viscosity (Poinsot 1989, Poinsot et al. 1991). The 
chemistry is represented through a single-step reaction A + B — ► P or a two-step 
process A + B — ► J, A + 1 — * P. The parameters of the Arrhenius- type reaction rate 
model chosen for the single-step simulations are close to that of global Methane-Air 
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combustion. The second mechanism is an attempt to study the effect of radical-like 
species in turbulent non-premixed flames. The second step (A + I — ► P) proceeds 
with an activation energy four times smaller and an enthalpy of reaction five times 
larger compared to the first step. 

The turbulent simulations are initialized by a well characterized one-dimensional, 
planar, laminar non-premixed flame which has been previously generated by the 
code (Chen et al 1992). A prescribed turbulent field interacting with the initially 
laminar flame is then allowed to decay. The value of the initial Taylor Reynolds 
number is the order of fifty, a value relevant for some jet flames (Dibble et al. 
1986; Magre et al. 1988), where the Reynolds number based on the velocity profile 
half- width increases from fifty to eighty along the jet axis. 

The Damkohler number, defined as the ratio of the turbulent eddy turnover time 
to a chemical time based on the heat release of the initial laminar flame, is varied 
over a range corresponding to fast and slow chemistry in order to study the effect 
of finite-rate chemistry on flame structure. The flame thickness is such that the 
flame-flow interactions are limited to a range of reactive and turbulent scales of the 
same order of magnitude. 

2.1 One-step chemistry model 

It is observed in the case of intermediate value of the Damkohler number (regime 
ranging between slow and fast chemistry) that points exist which are influenced 
by the fully burning and extinguished locations. In the plane (Ya,Z) defined by 
species A mass fraction and the mixture fraction (a conserved scalar, Z), these 
points are located between the lower bound on the species concentration distribution 
(equilibrium line, stoichiometric value, Zst = 0.5) and the upper bound (pure 
mixing), Figure 1. This is representative of a turbulence-induced mixing influencing 
both full burning and frozen flow points, and it corresponds to a transient response. 
Indeed, in the (Ya,Z) space, the concentration of a fluid particle is governed by 
(Lagrangian frame) : 


dY A j 7 

** £ 0®fS) 

According to this equation, depending on the time scale involved in the small 
scale diffusion process (enhanced by the turbulence) and in the reactive phenomenon 
(u^), the concentration of the fluid particle in the plane (F^, Z ) could reach all the 
values in the domain that lies between the pure mixing line and the fast chemistry 
line. The transient response is due to the effect of diffusion at small scale which is 
affected by the turbulence modifying the instantaneous species profile and, there- 
fore, the reaction rate. This transient response is representative of the finite-rate 
chemistry effect. 

It has been verified (Chen et al 1992) that in the databases, the stoichiometric 
surface and the location of peak reaction zone coincide, and also that contours of 
Z are parallel within the reaction zone. The mixture fraction dissipation rate \ 
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FIGURE 1. Distribution of the reactive species Y A versus mixture fraction Z. 

(X = 2D | VZ | 2 ), is found to be well correlated with the tangential strain rate 
evaluated at the flame surface. Gibson’s theory for constant density flows (Gibson 
1968) predicts a value of -0.5 for the correlation coefficient between strain rate and 
scalar dissipation rate. Incompressible DNS (Nomural et a l. 1992) report values 
ranging from -0.4 to -0.5. The present simulations provide a value of -0.6. This 
larger value is likely to be a result of dilatation associated with heat release in our 
computations. As expected, the scalar gradient is also found to be most probably 
aligned with the most compressive strain rate direction. 

In accord with laminar flamelet theory (Peters 1986), the scalar dissipation rate 
increases with reaction rate until a critical value is reached at which extinction 
occurs. Early in the simulation (i.e. within one eddy turnover time), the maximum 
value of the reaction rate interpolated along the local flame surface normal vector 
and plotted versus the scalar dissipation rate follows the common laminar flamelet 
response, Figure 2. However, when the flame is undergoing “full interaction” with 
turbulence, see Figure 3, a deviation from the bounds indicated by laminar flamelet 
is observed and is related to the transient effects observed in Figure 1. 

The reason for this deviation is due to the reaction rate being influenced by the 
local temperature but also by the species mass fraction. It appears that turbulent- 
enhanced mixing convects more species to the reaction zone than in a pure strained 
laminar flamelet. The vorticity does not penetrate the flame except at locations 
where extinction is observed. Indeed, within the flame region, the dynamic viscos- 
ity increases with temperature, creating a zone of large dissipation which has the 
effect of damping the turbulence. These observations suggest that even if the fluid 
mechanics features of the flow are close to the flamelet regime the dynamic infor- 
mation carried by the turbulence introduces some transient effects that certainly 
need to be included in the modeling to predict the mean burning rate capturing 
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1 / (X / Xlam) 

FIGURE 2 AND 3. Distribution of the maximum reaction rate with respect to the 
inverse scalar dissipation rate 1/ (x/x/am), for two successive times. 

with accuracy finite-rate chemistry effects. 

The pdf modeling shifts the problem of the evaluation of the mean burning rate 
to the estimation of the small scale mixing (Dopazo 1992). Hence, this approach 
is promising when the mixing process is not too much influenced by the reaction, 
allowing the chemical kinetics to decouple from the fluid dynamics. The distribution 
of the reactive species dissipation rate xy a (\y a = ZDa\ | 2 ) is a quantity 
representative of the intensity of mixing and is found in the simulations to be 
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FIGURE 4. Distribution of the reactive species A dissipation rate with respect to 
the mixture fraction Z , slow chemistry. 



FIGURE 5. Distribution of the reactive species A dissipation rate with respect to 
the mixture fraction Z, fast chemistry. 

sensitive to the Damkohler number. Actually, in the limit of very fast chemistry, 
the properties of the small scales of the reactives species field are known to be more 
sensitive to the chemical reaction than to the turbulent cascade energy process 
(Kuznetsov et al. 1990, p61). Thus it is hardly surprising that the reactive effect 
can still be perceived for intermediate values of the Damkohler number. 
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FIGURE 6. Distribution of the mixture fraction dissipation rate \ versus mixture 
fraction, slow chemistry. 

In Figure 4, two zones are observed, corresponding to Figure 1 or Figure 3. 
When the Damkohler number is large, Figure 5, the dissipation rate is negligible 
for Z < 0.4, whereas for the slower chemistry, Figure 4, a finite bound appears. 
This figure displays the penetration of the species A in the other side of the flame, 
creating, at the very least, a partially premixed situation. Here, the same physical 
turbulent-induced processes are observed than those previously noticed in the inves- 
tigation of the local flame topology. Locally, the gradient of the reactive species are 
simultaneously influenced by both reaction rate and turbulent transports. There- 
fore, even in non-premixed combustion, it is difficult to conclude that small scale 
mixing is unaffected by the reactive processes. 

Consequently, the estimation of a mean mixing time scale (used in any mixing 
model) without including information related to the reactive activity can lead to 
large disagreement between numerical predictions and measurements since the fluid 
dynamics motion is then somewhat disconnected from the kinetics (Chen et al 1990, 
Borghi et al 1990). A way to overcome this deficiency could be the introduction in 
the modeling of informations related to the gradient of the species field (Fox 1992a, 
Mantel et al 1990). 

The same approach is used to study the sensitivity of the mixture fraction dis- 
sipation rate x to finite-rate chemistry and turbulence-induced transient effects 
connected with those mentioned above. A detailed comparison of reactive and non- 
reactive fields indicates that the distribution of the mixture fraction dissipation rate 
is sensitive to the chemical processes. This is attributed directly to the modification 
in density and dynamic viscosity as a result of heat release. 

The distribution of the mixture fraction dissipation rate \ appears to lie within 
an envelope that at first sight follows the response of a typical laminar-like flamelet, 
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Figure 6. This has also been observed by Fox (Fox 1992b) in the study of the 
layer-like lamellar system. Usually, it is assumed for modeling purposes that the 
mixture fraction and its dissipation rate are uncorrelated quantities (Kuznetsov 
1990 p60, Peters 1986, Warnartz et al. 1986). Figure 6 indicates that the turbulence 
randomizes the values but that they are not completely uncorrelated. 

2.2 Two-step chemistry model 

The structure of the reaction zone obtained with the two-step chemistry in the 
case of slow chemistry is more complex. The intermediate species, I , field is such 
that a production-recombination zone lies on the oxidizer side of the domain (Z < 
0.5), and a diffusion zone is observed on the fuel side, Figure 7. If the second 
reaction was suppressed, the dissipation rate of the intermediate species \Y, (xy, = 
2T>i\ VYj | 2 ) would present two symmetric peaks on either side of the value Z = 0.5 
with a minimum at this point. Since / is consumed by A in the second reaction 
step, the maximum value of the dissipation rate is enhanced relatively to the peak 
on the B side, with a minimum shifted to approximately Z = 0.4 (location of the 
generation of I by the first step), Figure 8. 

The two reaction zones are not completely segregated in physical and mixture 
fraction space. However, for the same value of the Damkohler number, the global 
contribution of the reaction to the energy source term is broader in mixture fraction 
space than the corresponding contribution in the case of one-step chemistry. This 
event makes the flame less susceptible to undergoing local extinction compared to 
a flame modeled with single-step. This suggests that the modeling of extinction is 
strongly connected to the choice of the chemical scheme in terms of number of steps 
and species involved. 

3. Unequal diffusivity effects in non-premixed flames 

Within hydrogen or hydrocarbon flames, mono-atomic hydrogen and diatomic 
hydrogen atoms are present with heavier species that diffuse more slowly. The 
laminar flame structure is known to be strongly dependent on the larger diffusivity 
of radical species, which plays an important role in the ignition processes. One way 
to determine if, despite the large turbulent-induced transports, unequal diffusivity 
effects persist in turbulent non-premixed flames is to compare, from experimental 
results, the value of mixture fraction based on hydrogen atom with those based on 
carbon atom. Vranos (Vranos et al. 1992) investigated Methane-Hydrogen flames at 
high Reynolds numbers and reported differences between the two different mixture 
fractions, indicating that unequal diffusivities effects are present. 

In most of the turbulent non-premixed combustion models, the hypothesis that 
all the species diffuse at the same rate is used. Thus the analysis of DNS databases 
including Schmidt number effects in non-reacting and in reacting flow is relevant. 

3.1 Unequal diffusivity effects in non-reacting flows 

To investigate the effect of varying the Schmidt number on the small scale mixing, 
first DNS of non-reacting flows have been carried out. 


418 L. Vervisch 


0.25 



0 0.2 0.4 0.6 0.8 1.0 

Mixture Fraction 


FIGURE 7. Distribution of the intermediate reactive species I versus mixture 
fraction, slow chemistry. 
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FIGURE 8. Distribution of the dissipation rate of the intermediate reactive species 
I versus mixture fraction, slow chemistry. 

In the free decaying turbulence configuration, a species C is added on the “A- 
side” with a Schmidt number equal to one-half. The laminar initial pure diffusive 
zone is then allowed to interact with the turbulence. If the scatter plot of Yc versus 
Ya does not show the evidence of differential diffusion effect as in the experiment of 
Smith (Smith et al. 1992), the distribution of the ratio r\ = Figure 9, clearly 
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FIGURE 9. Distribution of the ratio r* = versus Ya\ pure mixing situation. 



Figure 10. 


Ratio r 2 = 


1 VV ^I J Ya versus Ya\ pure mixing situation. 
\vy c \ 7 \y a =*a) 


demonstrates the presence of this differential diffusion effects. Figure 9 displays for 
the same time the laminar equivalent evolution (where the exact solution is known) 
and the scatter plots of r x corresponding to the turbulent flow (the equal diffusivity 
problem would provide r x = 1). The turbulence induces a spreading of the points 
around the laminar response, and when Ya tends to zero, the ratio r i also tends to 
zero. Both of these observations reveal that small scale diffusion is sensitive to the 
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FIGURE 11. pdf of the differential diffusion variable; pure mixing situation. 


smaller Schmidt number of the species C. The ratio of the conditional mean values 
for the magnitude of the gradient of the species A and C, r 2 = 1 is 

presented in Figure 10. This quantity reaches its maximum when Ya — 0.5 and its 
minimum at the edge of the composition space, with large variations around the 
value r 2 = 1 corresponding to equal diffusivity. Kerstein (Kerstein 1990) proposed 
to define a differential diffusion variable, /, as: / = Ya — Yc- The pdf of this 
quantity is proposed in Figure 11. This pdf is somewhat symmetric around / = 0, 
with unequal zero value tails. The examination of the three previous quantities 
demonstrates that differential diffusion effects occur in pure mixing simulations. 
The turbulence induces a deviation from the laminar diffusive motion, but the 
unequal diffusive effect can still be perceived. 

Increasing the Reynolds number will decrease the scales at which the laminar 
diffusive process is acting. In large Reynolds number reacting flows, because of 
significant increase in kinematic viscosity with temperature, the differential diffusion 
effects pointed out above will certainly still be present at the diffusive scales where 
chemical activities occur, especially for very fast diffusive radical species. 


3.2 Unequal diffusivity effects in reacting flows 

If the Schmidt number of all the species is assigned a value of unity, the pdf 
of the curvature of the three dimensional flame is found to be symmetric (the 
surface exhibits both positive and negative curvature). However, when the Schmidt 
number of the species A is changed to one half (one-step chemistry) the probability 
of negative curvature along the stoichiometric surface increases, Figure 12. Negative 
curvature in this instance corresponds to reaction zones that are curved into the 
“B-side”’. This observation is correlated with the broader reaction zone found on 
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this side in mixture fraction and physical space. The same trend has been found by 
Trouve (Trouve 1991) in the case of premixed combustion. 



FIGURE 12. Pdf of the curvature along the flame surface; fast chemistry. 

In the two-step chemistry case without changing the initial turbulent field but 
modifying the Schmidt number of the intermediate species, I, DNS allows us to 
study the impact on the flame of a fast diffusive radical-like species. The contour 
plots of the intermediate species dissipation rate presents two areas related to the 
production-consumption and pure diffusion zones of I. For Scy t = 0.5, the two 
zones are wider than for Scy, = 1, Figure 13. This observation is related to the 
penetration of the rapidly diffusing intermediate species into both sides of the flame, 
which modified the magnitude of the dissipation rate but also the shape of the joint 
pdf of the dissipation rate and mixture fraction (compare Figure 8 and 14). 

These results suggest that globally small scale mixing is sensitive to the fast 
diffusive intermediate species, even though turbulent diffusion is the more active 
process in term of mean statistical description of the reactive flow field. Further 
examination of the non-premixed flame response to the unequal diffusivity effects 
are in progress. 

4. Modeling of non-premixed flames 

4-1 Introduction 

Usually, in industrial computer codes, transport equations for mean quantities 
(mean mass fraction, mean velocity and mean energy) are solved with the help of 
classical turbulence models for the Reynolds stress and the turbulent flux (usually 
Favre averaging is chosen, A = )• The validity of these turbulence models in 

the case of reacting flows with variable density will not be discussed here. The 
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FIGURE 13. Intermediate species dissipation rate contours, on the left Scj = 1, 
on the right Scj — 0.5. 
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Mixture Fraction 

FIGURE 14. Distribution of dissipation rate of intermediate species I versus 
mixture fraction; Scj = 0.5. 

mean reaction rate that appears in the mean species transport equation needs to 
be modeled and the aim of the method proposed in the following is to built a 
kind of “pdf generator” which generates, at a low cost, an approximated pdf of 
the thermochemical variables. This is with enough accuracy to compute in non- 
premixed flame in air the mean reaction rate of the oxidizer, woi* 

The chemistry is described through one global reaction step and heat losses by 
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radiative transfer are neglected; therefore, two variables (mixture fraction and oxi 
dizer mass fraction) only axe required for the description of the reactive flow field. 
This is a very crude approximation, but for engineering purposes, instead of using 
the local equilibrium hypothesis, it is useful to compute a first approximation of 
the reactive flow field which accounts for the interaction between the turbulent and 
chemical processes. This approximation can then be used to study the behavior of 
global parameters or to initialize more precise and more time consuming computer 
methods, such as Monte Carlo Methods able to deal with complex chemistry (Jones 
et al. 1985). 

The DNS results of turbulent non-premixed flames discussed previously indicate 
that transient effects have to be included in the modeling. Especially, the scatter 
plots of the reactive scalar concentration, temperature, and reaction rate show that 
the pdf of the reactive species exhibits a continuous shape between the frozen flow 
line and the equilibrium line. 

Two curves yb x = ybz( Z ) an ^ Vox = Voz( Z ) clearly limit the spreading of the 

pdf in the plane (Y 0x ,Z), {y l Qz < Vox V z )- The mean burnin 6 rate can be 
computed from an estimated pdf P(Yoz,Z ) supported by an area located between 
VOz an< l ybz • Three pieces of information are needed to perform this evaluation, the 
approximated location of the pdf in the concentration space, the extend of spreading, 
and the local amplitude of the pdf. The location of the pdf and the control on the 
extend of spreading will be based on the simplified probability density function 
approach (Borghi et al. 1983, Gonzalez et al. 1986). A way to evaluate the local 
amplitude of the pdf is derived from the DNS observations. 

4-2 A new formulation for the presumed pdf method 

Using information issued from the mean flow field, Borghi (Borghi et al. 1983) 
first proposed to built an approximation of P(Yoz, Z). This approximated pdf was 
obtained with the help of the Dopazo (Dopazo et al. 1976) small scale mixing 
model (also call iem Lagrangian model) and possessed some of the characteristics 
of the mean reactive scalar field. This model included some properties of the con- 
ditional pdf T c (Yoz\Z), (T c {Yoz\Z) = P %°zf } ), where P( z ) was Presumed as a 

Beta function with the two parameters Z and Z' 2 . 

Inside a two-inlet homogeneous reactor, when the pdf is discretized with two 
particles (one for each inlet), the evolution of the particles concentration equivalent 
to the pdf transport equation closed with the iem mixing model is (Pope 1979): 


dyoz _ d>o z Ti 0x + Yqx — yox 
dZ ~ Z-Z 

The quantity r t0z is a mixing time scale, and when only mean values are known, 
the previous equation is a closed form of Eq.(l). According to DNS observations, 
a transport equation must be solved for the mean dissipation rate of the re active 

< 2 y' J 

species £y 0i and also for the rms of the reactive species Y Qx to get r t0i = 
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This estimation of r tox allows the inclusion of information coming from the reactive 
processes into the mixing model. 

The homogeneous situation related pdf, P c (Yo x \Z), can be written: 

T c (Yo'\Z) = 6(y i S?(Z)-Y ox ), (3) 

where y l S™{Z) is composed of two trajectories obtained by solving (2) with the 
boundary conditions: {yox = 0, Z = 1) and (yoz = 1,Z = 0). 

In non homogeneous flows, the shape of the conditional pdf, P c (Yox\Z ), spreads 
out, and Borghi (Borghi et al 1983) proposed a way to approximate this pdf 
in premixed flames. This method has also been extended to non-premixed flame 
(Vervisch 1992). The pdf, presumed with the help of these previous formulations, 
reproduced only the mean value of the reactive flow field and was obtained in a 
discontinuous form, including large holes in the concentration space; moreover, the 
amplitude of the reactive species pdf was deduced from the amplitude of the mixture 
fraction pdf. These features are in disagreement with measurements (Masri 1988, 
Cheng 1991) and previous observations from DNS. An attempt to improve these 
presumed pdf models is now proposed. 

In the DNS field, if the conditional pdf P c (Yox\Z) effectively depends on Z, the 
pdf P c (£ |Z) where £ = \y~Z V J? £ ls only weakly dependent on Z, (Figure 15). 

*Ox VOx 



FIGURE 15. Conditional pdf P C (£|Z). 

Therefore, in the plane ( Yo x , Z ), each line between y}y x and y\> z% for a fixed value 
of Z, can be weighted with a reference pdf, /(£). The function /, able to mimic 
situations ranging from low spreading (low oxidizer fluctuations) to the bimodal 
form (large oxidizer fluctuations), has to be defined as: 
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Pc(Yox\Z) = 0 if Yox <ybz( z ) or Yox > ybx( z ) 

P c (Yox\Z)dY ox = f(Od( (4) 

f ° x(Z) P c (Yox\Z) dYox = f f(0 d£ = 1 VZ . 

Jy'o'(Z) Jo 

The estimation of the position of the pdf in the oxidizer mixture fraction plane 
is achieved when the curve y x Qz (Z ) and y 2 0z (Z) are known. The iem trajectory 
(which contains information on the join t effe ct of small-scale mixing and reaction) 
is taken as a reference. The quantity yg- = /„ y'6?( z )P( Z ) dZ is computed^ 
The following possibilities may be considered to determine the curves y Q z (Z) and 
ybx ( z ) (B° r shi et al. 1983): 



vbx = y l ox( z ) 

vl, = &.W 


Equilibrium Line (r <0i — > oo) 

iem Line 

y'6? < Y ox 

iem Line 

Mixing Line (r t0l . — > 0) 


The estimated boundaries being determined, if /(£) is a normalized Beta function, 

f(fi\ = . then the constraints, 

a' 

Jo 


Vox = f [ YoxP(Yox,Z ) dYox dZ 
J 0 Jo 


ylJZ) _ 

YoxPc(Yox\ z )dYox]P(Z)dZ 

v'oJZ) 


') 


-u£ 

Y' 0z = f [\y 0 x -Yox) 2 P(Yox,Z) dYox dZ 

Jo Jo 


(5) 


-as: 


vbJZ) 


ybjz) 


(Yox -Vox) 2 Pc(Yox\Z) dYox P ( Z ) dZ 


) 


according to the properties of the ft function, lead to: 


(gi -l)(l-ft2) 

flj(Q 2 — 1) 

= (l-fta)fii -n 2 (l-fti) ’ 


( 6 ) 


where : 
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ft , = r°,-^o, 

y 2 ox - y l o z 

^Oi — 2 ( ybziybx ~ vbx) ~ ( y 2 ox — i^ox)^ o*) ^1 
{vox -ybz) 2 Vi 
((y^) 2 ~2 ¥o*Y 0z + Y 2 0t ) 

(vbx -ybx ) 2 ^ i 


(In the equations, the mean values are defined as: a = a(Z)P(Z) dZ ,) 

The mean reaction rate &>ox can then be computed to advance the equation in 
time in the Navier Stokes solver, 


wox = f f uoxP{Yot,Z) dYoz dZ 
J o Jo 

■id: 


vlAZ) 


y'o z (Z) 


uoxPc(Y 0x \Z) dYox P(Z ) dZ 


) 


( 8 ) 


this is performed until the steady state is reached. 

This dynamic “pdf generator” is able to follow different shapes of pdf according to 

the level of turbulence (through the quantities, Z, Z' , Y Ox , Yox 5 T t 0 *i chemistry 
model parameters). The finite rate chemistry effect is implicitly included in the 
mixing model, especially when the time scale includes informations related to the 
chemistry. Unlike the previous formulation, the pdf is presumed in a continuous 
form, reproducing the transient effect observed in the DNS simulations. This simple 
model can be improved using more elaborate mixing model (Gao 1992) and can be 
of great interest to initialize the Monte Carlo method to reduce the large amount 
of computer time. 


5. Future plans 

The study of the turbulent flame structure modeled with two-step chemistry will 
be completed introducing a third step, B + 1 —► P, allowing the species / to undergo 
recombination on both sides of the flame. 

The configuration also needs some improvement; a shear will be included in the 
simulations using the methodology developed by Trouve (Trouve 1992). This is to 
study a model flow problem with properties closer to the non-premixed jet flame 
situation. 
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